Astronomy & Astrophysics manuscript no. 15423tex 


©ESO 2010 


September 8, 2010 





Letter to the Editor 



Dynamics and stability of relativistic GRB blast waves 

Z. Meliani 1 , R. Keppens 1 ' 2 ' 3 

1 Centre for Plasma Astrophysics, Celestijnenlaan 200B, 3001 Heverlee, K.U.Leuven, Belgium 

2 FOM-Institute for Plasma Physics Rijnhuizen, Nieuwegein, The Netherlands 

3 Astronomical Institute, Utrecht University, The Netherlands 

e-mail: zakaria . meliani@wis . kuleuven . be , rony . keppens@wis . kuleuven . be 

Received ; accepted 

ABSTRACT 

Aims. In gamma-ray-bursts (GRB), ultra-relativistic blast waves are ejected into the circumburst medium. We analyse in unprece- 
dented detail the deceleration of a self-similar Blandford-McKee blast wave from a Lorentz factor 25 to the nonrelativistic Sedov 
phase. Our goal is to determine the stability properties of its frontal shock. 

Methods. We carried out a grid-adaptive relativistic 2D hydro- simulation at extreme resolving power, following the GRB jet during 
the entire afterglow phase. We investigate the effect of the finite initial jet opening angle on the deceleration of the blast wave, and 
identify the growth of various instabilities throughout the coasting shock front. 

Results. We find that during the relativistic phase, the blast wave is subject to pressure-ram pressure instabilities that ripple and 
fragment the frontal shock. These instabilities manifest themselves in the ultra-relativistic phase alone, remain in full agreement with 
causality arguments, and decay slowly to finally disappear in the near-Newtonian phase as the shell Lorentz factor drops below 3. 
From then on, the compression rate decreases to levels predicted to be stable by a linear analysis of the Sedov phase. Our simulations 
confirm previous findings that the shell also spreads laterally because a rarefaction wave slowly propagates to the jet axis, inducing 
a clear shell deformation from its initial spherical shape. The blast front becomes meridionally stratified, with decreasing speed from 
axis to jet edge. In the wings of the jetted flow, Kelvin-Helmholtz instabilities occur, which are of negligible importance from the 
energetic viewpoint. 

Conclusions. Relativistic blast waves are subject to hydrodynamical instabilities that can significantly affect their deceleration prop- 
erties. Future work will quantify their effect on the afterglow light curves. 
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1. Progress in GRB dynamics modeling 

Gamma ray burst (GRB) outflows represent the most extreme 
variant of relativistic jets. That they represent finite opening an- 
gle flows is supported by the appearance of jet breaks in GRB 
afterglows ( Fra il et a l.||7997 ) (insights provided by a coupled 
simulation and light-curve synthesis point out subtle chromatic 
effects, van Eerten et al. 2010). Other convincing arguments rely 
on the polarization in prompt emission (Lazzat i et al.|2004| ), as 
well as theoretic al reasoning to invo ke physically plausible ex- 
plosion energies (Bloo met al.|2003| ). 

However, despite the progress in GRB understanding al- 
lowed by current detailed observations, various questions re- 
main. In particular, for the (micro and macro) physics occur- 
ring within the relativistic shock fronts. To improve constraints 
on GRB progenitor stars and overall GRB energies, we need to 
develop improved modeling, targeting the interaction of colli- 
mated relativistic blast waves with circumburst progenitor sur- 
roundings. This interaction is subject to non-linear physics af- 
fecting the details of how energy is transferred from the colli- 
mated ejecta to the circumburst medium. Studies of this inter- 
action rely on accurate, shock-capturing high resolution simu- 
lations. Previou s ID studies using Lorentz factor 100 shells in 
either uniform Meliani et al. 2007 ) or fully wind- structured sur- 
roundings (Meliani & Keppens 2007 ) quantified the need to re- 
solve a scale ratio of up to 6 orders of magnitude between shell 



width and traversed distance, to study the blastwave until the 
Sedov phases. 

Pioneering work in GRB dynamics in afterglow phases 
was perfomed anal ytically, e.g. |Meszaros & Rees (1997]), and 
numerically, e.g. Kobayashi et al. " fl999| ), using isotropic ex- 
plosion models. Rhoads (1999]) studied GRB outflow in two- 
dimensional settings analytically, by making assumptions anout 
sideways jet-expansion speeds. In the past decade, two and even 
three-dimensional relativistic simulations emerged, necessarily 
sacrifi cing in resolution (Gran ot et al. 200 1| Cannizzo et al.| 
2004 |, grid-adaptive studies at very high effective re solution 
( [Zhang & MacFadyen| [20091 |van Eerten et al.||2010] ) as well 
as high initial Lorentz factors (up to y = 100 Meliani et al. 
2007) became feasible. In contrast to earlier analytic assump- 
tions, these simulations found a slow lateral expansion of the 
GRB jet, in accord with semi-analytical works ([Kumar & Granot 



2003 ). However, the stability properties of the decelerating rel- 
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ativistic blast wave have not yet been addressed, in contrast, 
to the Newtonian regime, where the self- similar Sedov blast- 
wave structure is known to be unstable ( Vishniac 1983 ) when 
the compression rate is high (in turn a function of polytropic in- 
dex r < 1.3), the relativistic regime is more complex since the 
compression rate depends on effective polytropic index, as well 
as on the Lorentz factor. Relativistic blast-wave stability has to 
be studied in a single dimension, and can only be performed at 
extreme resolutions, as demonstrated here. 
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Fig. 1. The density variation at local time t = 66.9 days, with the large-scale view in the background, and a zoom into the shell at the 
bottom right (zooming into the square at the bottom of the large-scale view), yet further into zoomed to the shell edge at top right 
and in bottom left (zoomed on the squares in the full shell view). In the zoom figures we use inverse colors. 



2. Model setup 

The two-dimensional hydrodynamics blast wave simula- 
tion in this paper was performed with the special rela- 
tivistic hydro-module of the adaptive-mesh refinement code 
AMRVAC ( [Meliani et aLl|2007|). W e use a Synge-type equa- 
tion of state (Meliani et al. 2004) that takes full account of 
the change in effective polytropic index. In the present paper 
we neglect the effect of magnetic fields and ra diative losses. 
As an initial condition, a self-similar relativistic Blandford & 
McKee| (1976]) blast wave structure, given by their equations 
Eqs (28)-(30), is assumed to have a half-opening angle 6 = 20°, 
energy £j et = E{ so 6 2 /2ergs, and equivalent isotropic energy 
E{ S0 = 2.6 x 10 51 ergs. The blast wave propagates in a constant 
density medium n = 0.78cm" 3 . We begin the simulation at a 
shock Lorentz factor y s h = 25 (implying a local Lorentz factor 
of y 25/ V2 - 17.7), making the initial radius of the blast 
wave Rq ^ 2.68 x 10" 2 pc at local time ^ OC ai,o - 51.68 days 
and observer time £ bs,o = Ro/2y 2 c ^ lhour. At that time, 
the blast wave energy is mostly concentrated in a thin layer 
AR ^ R ly 2 = 8.5 x 10" 5 pc. To follow the evolution of this blast 
wave from relativistic to Newtonian phases, we use a spherical 
domain spanning [1.57 x 10" 2 , 1.57] pc in radius and [0, 45°] in 
meridional angle. To simulate this blast wave over this large dy- 
namical range, we need adaptive mesh refinement (AMR). This 
AMR run exploits 13 grid levels and uses an error estimator 
for normalized second derivatives on density with a weight of 
0.9 and of pressure with a weight of 0.1 as criterion for refine- 
ment/coarsening. In addition, we enforce coarsening in the back 
of the shell, and begin with 280 x 28 grid cells at level 1. By 
using a refinement ratio 2 between levels, we achieve an effec- 
tive resolution of (1.14688 x 10 6 )x (1.14688 x 10 5 ). Our block- 
tree AMR automatically tracks the shock front, while the num- 
ber of blocks remains nearly constant during the run: we have 
typically 8 blocks at the lowest level and about 50 000 blocks 
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Fig. 2. The variation in the Lorentz factor at the front shock 
against angle 6, at various times in the blast wave evolution. 

at the highest level. With this extreme resolving power, we can 
capture radial details at SR ^ 1.37 x 10" 6 pc. In the meridional 
direction, we have 66 - 0.0004°, and the initial shell contains 
50972 meridional cells. Figure [TJillustrates all the scales that en- 
ter this simulation. We used a predictor-corrector scheme with 
the HLLE solver (Harten et al. 1983) and minmod limited lin- 
ear reconstruction. This scheme avoids potential numerical in- 
stabilities at strong shocks (Mulle r et al.|2010|), as co nfirmed by 
analysis of its linearized form (Pandolfi & D' Ambrosio 2001| ). 

3. Lateral spreading and global deceleration 

During a long phase, the outflow propagates mainly spherically 
because of the relativistic initial radial speeds. However, the 
pressure contrast between the hot shell and the external medium 
causes a lateral spreading at the edge. In the meridional direction 
this consists of (1) a shock propagating in the external medium 
and (2) a rarefaction wave propagating in the shell towards the 
axis. Since the BM shell has a radially stratified pressure that 
declines sharply behind the front shock, the lateral flow expands 
mainly from the region near the head of the shell. The spread- 
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variation in the radial position of the frontal shock, 
to show a maximum unit value at each time, as a 
angle from axis. This shows the front deformation 
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Fig. 4. The variation in meridionafspeed at the shock front with 
angle 0, al several times. , , . . 

mg speed decreases quickly from just behind this region, in- 
ducing Kelvin-Helmholtz instability in the wings of the shell. 
This is fully consistent with the results by Zhang & MacFadyen 
(2009), and a zoom into these details is provided in Fig. [1] This 
Kelvin-Helmholtz development does not affect shell dynamics, 
since the characteristic time for the lateral flow is longer than the 
timescale found for the radial deceleration. 

The rarefaction wave that forms propagates inside the shell 
at the fluid's rest-frame sound speed. This wave extracts energy 
gradually, first from the edge and ultimately from the axial re- 
gion. During the relativistic phase, the shell becomes structured 
meridionally, and energy, pressure, and density ultimately de- 
crease from axis to edge. Deceleration occurs on the meridion- 
ally growing part of the shell affected by the wave, as quantified 
in Fig. [2] This meridional stratification of the blast wave makes 
the edge of the shell decelerate faster than the inner part, bending 
the shell away from its spherical shape, as shown in Fig. [3] This 
deformation becomes pronounced when the radial speed drops 
to about of the lateral speed, and increases the blast- wave work- 
ing surface and the distance traversed by spreading matter. This 
is all i n full agreement with previous simulations ([Granot et al. 
200Tj |Cannizzo et al.||2004| |Zhang & MacFadyen||2009|). The 



rarefaction front reaches the axis after t\ oca \ ^ 1 15.536 days. By 
this time, the blast- wave has decelerated to y = 1.09 at the axis 
(near-Newtonian phase). In Fig. [4] we quantify the meridional 
speed at the front of the shell for several times, to illustrate the 
propagation of the rarefaction wave. The meridional stratifica- 
tion of the blast wave decays when the Lorentz factor decreases 
to the point that causality links all parts of the shell together. 
In this near-Newtonian phase, energy becomes uniformly dis- 
tributed throughout the shell. 

4. Instability development in relativistic shocks 

The most novel discoveries of our simulation are from the de- 
tails of the initial evolutionary phase, where we now demonstrate 



that during the blast- wave deceleration, a ram pressure instabil- 
ity develops. This instability is found only while the blast wave 
speed remains relativistic, and disappears when its Lorentz fac- 
tor becomes y < 3 at time t\ 0C2L \ ^ 1 87 days (or observer time 
^obs - 53 hours, as determined from A^ b S = A^ oca i - AR/c). 
Figure [5] shows its spontaneous appearance and disappearance, 
manifesting itself during some 100 days (local time). This in- 
stability induces small-scale bending of the blast front at wave- 
lengths 86 ^ 0.2° (500 times larger than our cell meridional 
size, hence more than resolved properly). Detailed analysis in- 
dicates that each local front merely affects the minute region 
that is causally linked at this time, with causality angles around 
^causality -7° initially that increase during the evolution. The ra- 
dial size of these ripples is given by R86 ^ 10" 4 pc at local 
time 67.5 days. This is smaller than the thickness of the shell 
8R 3.7 x 10" 3 pc. The instabilities develop as small-scale 
ripples at the shock front, and the higher thermal pressure in- 
side them pushes the newly shocked external medium sideways 
to their left and right upstream. Fluid within the shell (behind 
the shock) is pushed downstream. This local structure generates 
vorticity in the post shock region, and is a prime candidate that 
can contribute to magnetic field amplification. 

The instabilities depend, as in the classical case ( Vishniac 
|1983| ), on internal pressure within the shell being thermal and 
isotropic, while external pressure is mainly ram pressure in the 
direction perpendicular to the shock front. This induces rip- 
ples, which grow in size as time proceeds, leading to shell frag- 
mentation. This fragmentation appears at a local time offset 
<^iocai ~ H5 hours after the start of our simulation, which is 
of the same order as the sound crossing time of the shell. At 
this time the shell has already decelerated to a Lorentz factor 
of order 10 and the local causality angle is reduced to about 7 
degrees. This instability appears initially at the outer edge (see 
Fig. [5] bottom), which is where the rarefaction wave first affects 
the shell. The edge thereby decelerates more quickly than the 
centre, is first to bend away from spherical shape and to reach 
the lower Lorentz factors. This is seen in the bottom panel of 
Fig. |5J while the second panel upwards shows that later the en- 
tire shell develops local instabilities (even before the rarefaction 
wave has affected its evolution). The size of the fragments in- 
crease because of the spherical expansion. Moreover, the density 
and Lorentz factor contrast between fragments and surround- 
ings increase and more energy become concentrated in these 
the fragments. They appear as independent shells propagating 
downstream, until the Lorentz factor of the entire blast wave 
becomes smaller then y < 3, or when they are reached by the 
inward propagating rarefaction wave. As the Lorentz factor of 
the shell gradually decreases, as well as the compression rate 
r = yrp, me soun d crossing time for the shell increases. The 
shell ultimately returns to a configuration that is stable against 
these pressure-ram pressure instabilities. In this near-Newtonian 
phase, the fragments diffuse and disapp ear, in a way tha t is con- 
sistent with linear analysis predictions ( Vishni ac|1983| ). During 
this phase, the blast wave propagates adiabatically with effec- 
tive poly tropic index of T e ff ^ 5/3 > 1.3. As the fragments 
are forced to move, they form tails in the rear of the BM shell. 
Our result is the first to show clearly that some of the as yet un- 
expected instability modes in BM solutions exist, decay slowly 
as the system evolves. In this case, small-scale structure appears 
when the Lorentz factor is 10, and disappears when it reaches 3. 

As the rarefaction wave propagates towards the axis, it not 
only forces the fluid to spread towards the edge, but damp- 
ens these pressure-ram pressure instabilities. In the relativistic 
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phase, the instabilities persist all the way to the axis. The rarefac- 
tion wave reaches the axis when the Lorentz factor is already low 
throughout, and the instabilities have almost completely disap- 
peared. After that, the shell follows the BM solution. However, 
since the rarefaction wave extracts energy, the energy density in 
the shell becomes lower than in the equivalent isotropic case and 
is not uniformly distributed. All these effects produce to a faster 
overall shell deceleration. At t\ oca \ = 1 166.27 days, along the axis 
the speed is 25 % lower than the value obtain ed from the ID ana- 
lytical solution ( [Blandford & McKee|1976] ), and 31% lower 20° 
away from axis. 

5. Conclusions 

We have carried out an extreme resolution simulation of prop- 
agating relativistic GRB jets in a constant density medium. We 
have concentrated on dynamical aspects previously missed. We 
have confirmed that the interface between the shell and the ex- 
ternal medium develops a rarefaction wave propagating into 
the shell, which induces slow lateral spreading. The region af- 
fected by the rarefaction wave loses energy and decelerates more 
rapidly. From edge to axis this deceleration becomes gradually, 
more efficient. The shell deforms, increasing the work surface 
and inducing an earlier shell deceleration. Our main result con- 
cerns the development of a relativistic pressure-ram pressure in- 
stability at the front of the blast wave. This instability induces 
differences in the deceleration properties during the relativistic 
phase. It is a relativistic analogue of the instability discussed in 
Vishniac (1983), but in relativistic disguise where the compres- 
sion rate is high even though the shock is adiabatic. The com- 
pression rate is then proportional to the Lorentz factor of the 
blast wave. Its ultimate fate is to decay slowly and disappear 
when the Lorentz factor decreases to y < 3. Finally, it should 
be noted that because of the strong beaming of the nonthermal 
emission from the blast wave, the instability that forms at the 
shock front might be expected to have an effect on the shape of 
the observed light curve, since at the front shock this instabil- 
ity can induce instantaneous variation in Lorentz factor where 

Tmax ~ 3y m j n . 
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Fig. 5. The shell edge at various times, increasing from bottom 
to top. This shows the small-scale instability development and 
evolution. 



